INTEGRABILITY OF A t - J MODEL WITH IMPURITIES 
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Abstract 

A t — J model for correlated electrons with impurities is proposed. The impuri- 
ties are introduced in such a way that integrability of the model in one dimension 
is not violated. The algebraic Bethe ansatz solution of the model is also given and 
it is shown that the Bethe states are highest weight states with respect to the su- 
persymmetry algebra gl(2\l). 
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1. Introduction 



The quantum inverse scattering method (QISM) has lead to many new results in the 
study of integrable and exactly solvable systems. Amongst these is the fact that the t — J 
model for correlated electrons is integrable in one dimension at the supersymmetric point 
J — 2t with the supersymmetry algebra given by the Lie superalgebra gl(2\l). This was 
made apparent in the works |TJ, 0] where it was shown that the Hamiltonian could be 
derived from a solution of the Yang-Baxter equation. Also, solutions of the model were 
found by means of the algebraic Bethe ansatz. 

One attractive aspect of the quantum inverse scattering method is that one is allowed 
to incorporate impurites into the system without violating integrability. In this context, 
several versions of the Heisenberg chain with impurities have been investigated [^, [|. |5| . 
For the specific case of the t — J model this idea was first adopted by Bares whereby 
the impurities were introduced into the model by way of inhomogeneities in the transfer 
matrix of the system. Another possibility was explored by Bediirftig et. al. with 
impurites given by changing the representation of the gl(2\l) generators at some lattice 
sites from the fundamental three dimensional representation to the one parameter family 
of typical four dimensional representations which were introduced in || to derive the 
supersymmetric U model. 

Here we wish to propose a third method for introducing integrable impurities into 
the t — J model. This is achieved by replacing some lattice sites with the dual space of 
the fundamental three dimensional representation. A significant point here is that only 
recently have new Bethe ansatz methods been proposed in order to solve such a system 
because of a lack of a suitable (unique) reference state. Rather one is forced to work with 
a subspace of reference states. This approach has been developed in the works of Abad 
and Rios || 10] and has already been adopted in to find a Bethe ansatz solution of 



the supersymmetric U model starting from a ferromagnetic space of states. 
The Hamiltonian of this t — J model with impurities reads 

H = + T ~hi,i+lQi ~ -7-Qihi^i+i (1) 



i=l i&l 



where 

hi,i+l 



~ ^2( c l*d+i,v + 4+i, CT Ci, CT )(l - rii _ ff )(l - n i+ i _ ff ) 

(J 

+2(Si.S i+ i - ^nin i+1 ) + n ; + n i+1 - 1, 

Qi = ^ v{4,a C * - 4<V) (! - n i -a) (1 - n-a) 

+2(Si.S-^n in )-n+l 
4 

and periodic boundary conditions are imposed. Above c\± 's are spin up or down annihila- 
tion (creation) operators, the Si's spin matrices, the rij's occupation numbers of electrons 
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at lattice site i. The Aj are arbitrary complex parameters and I is simply an index set 
with elements in the range 1, 2, L. We make the assumption that if i G / then i ± 1 ^ / 
since otherwise extra terms are needed in the Hamiltonian for integrability. The operators 
without site labels in the expression for Q { act on the impurity space coupled to the site i. 
Note however that the interactions involving the impurity sites are three site interactions 
involving the sites i and i + 1 as well as the impurity. The local space of states for an 
impurity site has the basis 

IT), II), III) 

in contrast to the local spaces for the other sites which have bases 

IT), ID, |o> 

as is the case for a pure t — J model. The reason for this choice is so the Hamiltonian 
conserves magnetization and particle number. Finally we mention that the first term in 
eq. ([!]) is the Hamiltonian for the pure t — J model. We can recover this model from 
eq. ([!]) by taking the limit Aj — > oo for each i £ I. 

In this paper we derive the Hamiltonian eq. ([!]) by means of the QISM which guar- 
rantees integrabilty. We will also find solutions to the model using the algebraic Bethe 
ansatz. Finally we also show that the Bethe states which are obtained by this procedure 
are in fact highest weight states with respect to the underlying supersymmetry algebra 

2. Derivation of the Hamiltonian 

Recall that the Lie superalgebra gl{m\n) has generators {Ej}™^ satisfying the com- 
mutation relations 

[E], E{] = b)E\ - (-l)(ra+0)(M+M)aj£;* (2) 
where the ^-grading on the indices is determined by 

[i] = for 1 < i < m, 

[i] = 1 for m < i < m + n. 

This induces a Z 2 -grading on the gl(m\n) generators through 

[E]\ = [z] + [j] (mod2). 

The vector module V has basis {i/}™^™ with action defined by 

E)v k = 8 k jV \ (3) 

Associated with this space there is a solution R{u) G End(V <8> V) of the Yang-Baxter 
equation 

R 12 (u - v)R 13 (u)R 23 (v) = R 23 (v)R 13 (u)R 12 (u - v) (4) 
on the space V <g> V <8> V which is given by 

R(u) = I®I-~ y %2e)®4(-l)W. (5) 

i,3 
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We remark that eq. (f|) is acting on a supersymmetric space so the multiplication of tensor 
products is governed by the relation 

{a®b){c®d) = (-l) [b][c] ac®bd (6) 

for homogeneous operators b, c. 

The solution given by eq. (|5]) allows us to construct a universal L-operator which reads 

L(u) =I®I--Y, e )® E K- l ) m - ( 7 ) 

This operator gives us a solution of the Yang-Baxter equation of the form 

R l2 {u - v)L 13 (u)L 23 (v) = L 23 (v)L 13 (u)R 12 (u - v) 

on the space V <S> V <E> gl(m\n) which follows from the commutation relations eq. (§). The 
dual representation to eq. (|J) acts on the module V* with basis {vi}™^ 1 and the action 
is given by 

E}v k = -(-lf + ^5lv r (8) 
By taking this representation in the expression eq. (0) we obtain the following i?-matrix 

R*(u) =I®I + -Ye)®e l A-l) [m (9) 

hi 

giving the solution 

R 12 (u - v)R* 13 (u)R* 23 (v) = R* 23 (v)R* 13 (u)R 12 (u - v) (10) 

on V ® V ® V*. 

We wish to construct an impurity model with generic quantum spaces represented by 
V and the impurity spaces by V*. To this end take some index set / = {pi,p 2 , ....,pi}, 1 < 
Pi < L and define 

L 
i=l 

where 

Wi = V if i£I, 

Wi = V®V* if iel. (11) 

In other words for each % G J we are coupling an impurity into the lattice which will be 
situated between the sites i and i + 1. 
Next we define the monodromy matrix 

T(u, {X}) = R 01 (u)R 02 (u)....R 0L (u) 

where we have 

Roi(u) = Roi(u) for i ^ I, 

Roi(u) = R w {u)R* w ,{u- \i) for iel. 
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Above, the indices %' and i" refer to the two spaces in Wi (cf. eq. (|TT|)) and the \ are 
arbitrary complex parameters. A consequence of eqs. (f|, [H|) is that the monodromy 
matrix satisfies the intertwining relation 

R 12 (u - v)T 13 {u)T 23 {v) = T 23 (v)T 13 (u)R 12 (u - v) (12) 

acting on the space V <E> V <S> W . The transfer matrix is defined by 

t(u) = tr cr T(u) (13) 

where the matrix a has entries 

a) = (-l)Mb^ 
from which the Hamiltonian is obtained through 

(14) 

In this derivation we used the property 

Qi{—2hij+x)Qi = Qi 

which follows from the fact that Q projects onto a one- dimensional space. This simplifies 
the calculations and is one of the reasons why this model is much simpler than other 
impurity chains. From eq. (|E|) we conclude by the usual argument that the transfer 
matrix provides a set of abelian symmetries for the model and hence the Hamiltonian 
is integrable. In the next section we will solve the model by the algebraic Bethe ansatz 
approach. The explicit form of the Hamiltonian eq. ([I]) is given by making the following 
identification between the basis elements of V, V* and the electronic states: 



^ = IT>, 


Vi = 


ID, 


- 2 = l!>, 


v 2 = 


IT), 


^ 3 = |o), 


l>3 = 


IU> 



H = -2 — \n(v L T(u)) 
du 



3. Algebraic Bethe ansatz solution 



By a suitable redefinition of the matrix elements, the solutions ([5], 0) may be written 
in terms of operators which satisfy the Yang-Baxter equations (|], [10]) without Zjzf-grading 
(see e.g. (12|]). These operators read 

R(u) = 5>«®eJ(-l) 



[*][?"] _ t e \ 
u 3 



R (u) 



+ -4 

u J 



1) 



h3 
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and hereafter we will use these forms. In the following we will also need the /^-matrices 



3 2 
r(u) = Y{-l)^4®e^- -e\®ei 



i,3=2 
3 

u 



r*(u) = ^(-lp 1 e:®e^ + -e}®e}(-l 

i,j=2 



which belong to a gl(l\l) invariant (6- vertex) system. From this matrices we define the 
monodromy matrices 

t(v,{u}) = r i(v - Ui)r Q2 {v - u 2 )...r m (v - u N ), 
t*(v,{\}) = r* 1 (v-\ i y 02 (v-\ 2 )...r* l (v-\ l ). 

First we construct the Yangian algebra which has elements {^j("u)}^t". Relations 
amongst these elements are governed by the constraint 

R 12 (u - v)Y 13 (u)Y 23 (v) = Y 23 (v)Y 13 {u)R 12 {u - v) (15) 

where 



Y{u) = Y J ^®y^ 



By comparison with eq. (|12D we see that the monodromy matrix provides a representation 
of this algebra acting on the module W by the mapping 

Ti(Y l j (u)) k l = (-l)(MM+[i]K+HW) T jfc( u ). ( 16 ) 
Moreover the transfer matrix is expressible in terms of this representation by 
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T(U 



J2(-l) [l]+m MYl(u)) 



i=i 



The phase factors present above are required since the Yangian algebra is defined with a 
non-graded .R-matrix. In the following we will omit the symbol 7r for ease of notation. 
For a given {a} = (ai, a 2 , ai), a« = 2, 3 we define the vector E W by 



L 



w 
i=i 



where 



w l = v 1 for i ^ I, 

w % = v 1 ® t> a . for i = pj G I. 

Now set X = span {v^}. It is important to observe that the space X is closed under 
the action of the elements YJ(u), i,j = 2,3 which generate a sub- Yangian. We may in 
fact write 

Y;(u)v^ = tf { ^(u,{\})v^ 
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which follows from the fact that the YHu) i,j = 2, 3 act trivially on the vector v 1 in the 



3 

sense 



y 2 2 ( n y = Yfi^v 1 = v 1 

Yi( u y = vi^v 1 = o 

Setting 

SW({u}) = Yfi( Ul )Yf*(u 2 )....Yf N (u N ), Pi = 2,3 
we look for a set of eigenstates of the transfer matrix of the form 

= J2 S^({u})v^F( pM (17) 

where the Ffp a \ are undetermined co-efficients. We appeal to the algebraic equations 
given by eq. (|T5|) to determine the constraints on the variables U{ needed to force eq. (|T7D 
to be an eigenstate. Although many relations occur as a result of eq. only the 
following are required: 

Yl(v)Yl(u) = a(u-v)Yf(u)Y 1 \v)-b(u-v)Tf t (v)Y 1 \u) (18) 
Y]'{v)Y?{u) = Y?{u)Yj\v)r%{v-u)-b{v-u)Y?\v)Y«{u) (19) 
a{v-u)Y?{v)Y?{u) = Yf{u)Yf{v)r^ a ,{v-u) (20) 



with a(u) = 1 — 2/u and b(u) = —2/u. All of the indices in eqs. (p"8|, 19) assume only the 
values 2 and 3. Using eq. ( |18D two types of terms arise when Yf is commuted through 
Yf. In the first type and Yf preserve their arguments and in the second type their 
arguments are exchanged. The first type of terms are called wanted terms because they 
will give a vector proportional to <3> J , and the second type are unwanted terms (u.t.). We 
find that 

JV 

Yf(y)& = a(v) L Y[ a(ui - v)& + u.t.. (21) 

i=i 

Similarly, for i = 2, 3 we have from eq. flT5| ) (no sum on i) 

Yt(v)& = ^'>({n})^(t;)^ } } (^{n})^^ Q} + u.t. 

= S^({u})t^} } (v, {\})v^F^ a} + u.t 

= S^'H{u})?$%(v, K X})v^F( p>a} + u.t. 

where 

The contribution to the eigenvalues of the transfer matrix is 

Yl(v)V + {-l)WY*{v)& = E(- 1 ) M+W ^K}(^ {«. + u.t. 

i=2 

(22) 
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At this point we need to perform a second- level, or nested Bethe ansatz procedure to 
diagonalize the matrix 

where we have used the fact that = unless [j] = [{[3, a}}. The above matrix 

is simply the transfer matrix for a g/(l|l) invariant system acting in the tensor product 
representation of N copies of the vector representation with inhomogeneities {u} and I 
copies of the dual representation with inhomogeneities {A}. 

To diagonalize this matrix we construct the Yangian generated by 

3 

y( u ) = Yl e )®yi( u ) 

i,3=2 

subject to the constraint 

n 2 (w - v)yi 3 (u)y 23 (v) = y23(v)yi 3 {u)r 12 {u - v). (23) 

From the above set of relations we will need the following 

y 2 2 (v)yl(u) = a(u-v)y 3 2 (u)y 2 2 (v)-b(u-v)y 3 2 (v)y 2 2 (u), (24) 
Vl(y)yl{u) = -a(u-v)yl(u)yl{v)-b{v-u)yl{v)yl(u), (25) 

yl(v)yl(u) = ~ a ( iu ~% 3 2 (u)y 3 2 (v). (26) 

a{v — u) 

Proceeding similarly as before, we look for eigenstates of the form 



with the vector w given by 



Using (Hill) ^ follows that 
with 



Ti(y)4> = Ai(v)(p + u.t. 



TV M I M 

A i( v ) = II a ( v ~ u i) \\ a {lk - v) - JJ a(v - Xj) Y\ a(7fe - v) 

i=l k=l j=l k=l 

The unwanted terms cancel provided the parameters 7^ satisfy the Bethe ansatz equations 
(BAE) 

N I 

Y[a( lk -u i )=Y[a^ k -X j ), k = l,2,...,M. (27) 

i=l 3=1 

Combining this result with eq. ([H]) we obtain for the eigenvalues of the transfer matrix 
eq. © 

TV 



A(v) = a(v) L Y[a( Ui - v) + Ai(«). {21 
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Cancellation of the unwanted terms in ( |2"T| , p2"D leads to a second set of BAE which are 

N / \ M 

a(u h ) L TT ^llhl = _ FT a ( 7fc _ Uh ) fc = i, 2 , iV. (29) 

We will not give the details proving the cancellation of the unwanted terms but remark 
that the calculation is analogous to that given in M for the pure t — J chain. 

Making a change of variable u — > iu + 1, 7 — > 27 + 2, A — > i\ + I the BAE read 

'^±*V = TT Mh -^ + 2 ' n Mh -^-' > /> = i,....,iv, (30) 

u h -ij f£ U/» - «i - 2i £jj Ufc - 7* + « 

7fc - u i ~* iTifc-^ -* 

In the absence of impurities (limit I — > 0) we recover the form of the BAE first derived by 
Sutherland [[0| and later by Sarkar [jTJj] for the usual t-J model. Adopting the string con- 
jecture, or more specifically assuming that the solutions Ui are real or appear as complex 
conjugate pairs and the Xj are real, we find string solutions 

<g = < + «(« + 1-2/3), a = l,2,...,iV n , P = 1,2,..., n, n = l,2,... 

and the 7^ are real. The number of n-strings N n satisfy the relation 

N = Y,nN n . 

n 

As was shown in the papers [I], |2| two other forms of the Bethe ansatz exist which are 
obtained by choosing a different grading for the indices of the g/(2|l) generators. Recall 
that the above calculations were performed with the choice 

[1] = [2] = 0, [3] = I. 

Choosing 

[1] = 1, [2] = [3] = 

yields the eigenvalue expression 

N M I 

A(v) = -a(-v) L Y[a(v-Ui) + ]Ja(v - -f k ) JJa(Aj - v) 

i=l k=l j=l 

N M 

+ Y[ a{v - Yl a{lk - v) 

i=l k=l 

subject to the BAE 

M 

a(-Ui) L = Ylai'jk - u^, i = l,2,...,N 
fc=i 

II , ~ 7 i = -U a ^-^U n V h = l,2,...,M. 
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In the limit I — > we recover Lai's [HJ (see also ) Alternatively, choosing 

[1] = [ 3 ] = 0, [2] = 1 

yields the eigenvalue expression 



N 



M 



A(v) = a(v) L \\a(ui - v) + Y\_ a ( v ~ Ik) \\a{^ 3 

i=l k=l 

N M 

- jj J a(ui -v)Y[ a ( v ~ Ik) 



k=l 



with the BAE 



M 



a[Uj 



= Y\ a ( u i~7k) i = 1,2, N, 

k=l 

N I 

Y[a{ui--f k ) = Yla(Xj-^k) fe = l,2,...,M. 



i=i 



Finally, from the definition of the Hamiltonian eq. (14) we see that the energies are 
given by 

E = -2^-\n (v L A(v)) 
dv v 

Using the eigenvalue expression eq. we obtain 



v=0 



E 



- 1 



where the Ui are solutions to the equations ( |30| , |3T| ). 

4. Highest weight property 

Next we wish to show that the eigenstates constructed in the previous section are in 
fact highest weight states with respect to the underlying supersymmetry algebra gl(2\l). 
The highest weight property of the Bethe states has been proved for many models, such as 



the Heisenberg chain [|T7| and its generalized version |l8fl , the Kondo model []19], the usual 
t — J model 0, the Hubbard chain p(| 21], |22] and its gl(2/2) extension [p3| |. However, 
as far as we are aware it has never been shown before in the case where a subspace of 
reference states has been used in the Bethe ansatz procedure. 
Let us begin by considering 

= £ E!S^({u})v^F( p>a} . 
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By means of the nesting procedure we know that the co-efficients F?p a | are such that we 
have the following identification of states 

S^({u})v^F( M = yl( ll )yl( l2 )....yl( lM )w 

for a suitable solution of the BAE. By comparing eqs. (f7,|23,16) it is possible to determine 
algebraic relations between the elements of the Yangian algebra and the supersymmetry 
algebra. For our purposes we need the following 

[El yl{u)Y p = -yl{u)1 + y!(ur,(-lp. (32) 

Noting that E%w = it is evident that we may write 

M 

Elyl( ll )....yl( riM )w = J2 x h X h 

h=l 

with 

x h = 2/2(71) yl(jh-i)yi(jh+i)--y%(jM)w 

and the Xh some yet to be determined co-efficients. To find Xh we write 

2/ 2 3 (7i)....*/ 2 3 (7mV = II ~" (7fe ~ 7 ^ 2 3 ( 7fe )X fe 



where we have used eq. (|26|). Now by using the relations (^,^,^) and looking only for 
those terms which give a vector proportional to Xi we find that 

h ~ 1 a (7 7 ) / 1 M N M 

Xk = n a ( h _ — r ( n ~ x ^ n ^ ~ ^ ~ n ^ ~ ^ n ^ ~ ^ 

j=i wj 7h) \j = i kj L h i= i k ^ h 

which vanishes because of eq. Q2"7|). Thus we see that 

Ej& = 0. 

Next we consider the action of E\ on <E> J '. Using eqs. (|7 UT5|Jl6|) we find the commutation 
relation 

[ElY?{u)] = 5%Y}{u)-Y?{u). (33) 
As before, since we can write the general expression 

h,p 

where 

and for any vector {w} we have 

{ w -} = ( Wlj W2i w h -i), {w£} = (w h+ i, , w N ). 
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To calculate zu,p we begin by writing 
& = S^\{ Uh })Y^(u h )S^H{ut})v^F( M 



h-l 



i=i 



where we have used the relation eq. (p0|) . Now applying eq. fl33l) and using the relations 
(|18|,|1^) to determine the terms which give a vector proportional to Zh,p we find that 

(JV N M N 

a(u h ) L JJ a(iij - u h ) - a(u h - u,-) a(7 fc - u h ) 
i^h i^h k=l 

which vanishes as a result of eq. (E^) . We then conclude that 

E\& = 

which completes the proof that the Bethe states are g/(2|l) highest weight states. We 
observe that this property can also be proved for the other two choices of gradings in a 
similar way. 

5. Conclusions 



In this paper we have introduced a new integrable version of the t-J model with 
impurities. The model was solved through an algebraic Bethe ansatz method and three 
different forms of the BAE were derived. A proof of the highest weight property of the 
Bethe vectors with respect to the gl(2\l) superalgebra was also presented. A possible 
application of these results would be an analysis of the structure of the ground state and 
some low lying excitations of the model in the thermodynamic limit. 
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